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We show how it is possible to more than double the on-site interaction energy of neutral atoms 
in optical potentials by the technique of radio-frequency (rf) dressing, while maintaining interwell 
dynamics. We calculate Bose-Hubbard parameters for rf dressed optical lattices and arrays of rf 
dressed dipole traps. We show that decreasing the distance between wells, by the interpolation of 
wells confining different itif states, increases the interaction energy more than decreasing the height 
of the classically forbidden region between existing wells. The schemes we propose have negligible 
Landau-Zener losses caused by atomic motion; this was a dominant effect in the first experimental 
demonstration of the modification of an optical potential by radio-frequency dressing. 

PACS numbers: 03.75.Lm, 34.50.Cx, 37.10.Jk, 64.70.Tg, 67.85.Hj 



The study of complex nonlinear quantum systems is a 
major area of research in the fields of atomic and con- 
densed matter physics. A subject of current interest is 
the effect of the nonlinear interaction term on the be- 
haviour of ultracold atoms confined in optical lattices. 
The dynamics of ultracold bosonic atoms in such a sys- 
tem has been shown to be described by the Bose-Hubbard 
Hamiltonian [l[ 

h = -j J2 ota+a-S^fa -1 ) (1) 

<i,j> i 

where i and j denote lattice sites of a homogeneous lat- 
tice. The ground state of the Bose-Hubbard Hamiltonian 
passes from superfluid to Mott insulator as the parame- 
ters controlling U and J are varied' this behaviour has 
been demonstrated experimentally |2|. 

The on-site interaction energy, the Hubbard U, is the 
dominant parameter characterising interactions between 
ultracold atoms in optical lattices. As such it plays a 
pivotal role in phase transitions Q and entanglement 
[g| . The ability to generate complex entangled states 
has drawn interest to these systems for the purposes 
of quantum computing Q and quantum simulation 
Typically, the magnitude of U controls the purity of the 
resulting many-particle state, or the speed at which this 
entanglement may be generated [6( . It is likely that there 
are a significant number of 'impurity' atoms present in 
the Mott insulator states currently being made [1, Q; 
for applications such as quantum computing and simu- 
lation it is desirable to reduce the number of these im- 
purities as much as possible. Furthermore, by increasing 
the on-site interaction relative to the tunnelling energy, 
it may be possible to push optical lattice systems into 
new regimes, for example where the on-site interaction 
energy is greater than the band gap 
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In this paper we study ways to increase the on-site 
interaction energy, U , while still maintaining slow in- 
terwell dynamics, using the technique of radio-frequency 
(rf ) dressing of optical potentials. We study both optical 
lattices and arrays of highly focused laser dipole spots, 
both seen aspromising systems for quantum information 
processing [5|, Q • We restrict our study to dressing the 
optical potential along a single direction of the lattice. 
The atoms in these potentials will have far lower Landau- 
Zener loss rates than in the recent experiment [101 ] . and 
low collisional loss rates. For a number of specific ex- 
perimental configurations we calculate the enhancement 
factor for U at a particular interwell tunnelling rate. 

Rf dressing has been used in a variety of experiments to 
dress magnetic potentials [ll|, EH • Dressing of optical 
potentials has received much less attention, although re- 
cently both theoretical [3] and experimental [l(| studies 
have been carried out. A dominant feature of the re- 
cent experiment [l(| was high loss rates, attributed to 
non-adiabatic Landau-Zener transitions We should 
mention other methods which may be adapted to increase 
on-site interactions in a lattice include tunable scatter- 
ing lengths and Feshbach resonances [HI, EH , and other 
methods proposed to give sub-half-wavelength structure 
to optical lattices, for example Raman processes [I?], EH ■ 

The phase transition from the superfluid state ceases to 
be adiabatic when the inverse timescale 1/r becomes of 
the same order of magnitude as the frequencies of the low- 
est lying excitations, which are ~ U in the limit U 3> J; 
increasing U will decrease the number fluctuations in the 
final state. Furthermore, when finite temperature effects 
are taken into account, it has been shown that increas- 
ing U increases the purity of the final state [l!| , as one 
would expect from thermodynamic considerations. Both 
temperature and non-adiabatic defects are exponentially 
suppressed by increasing U, as discussed below. 

We treat the case with the spin dependent potential 
only along one direction of the lattice; dressing a spin 
dependent potential of more than one dimension tends 
to lead to structures with extended potential minima 
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FIG. 1: (Colour online) The types of RF dressed potentials considered in this paper. Dressed potentials are solid 
black lines, undressed potentials are broken lines {nip = +1/ — 1 (green/red dashed), (blue dash-dotted)]. The 
Wannier functions for the lowest state are displayed below the potentials for two of the wells. The rap = +1 state of 
(a) is at an energy substantially higher than the states shown. The optical lattices of (a) and (b) have a + 
polarisation. The recoil energy is defined as E r = h 2 /2m\ 2 . Parameters: 
(a) A = 790.06 ran, B Q = 300 G, v rf = 213.65 MHz, Q rf /2Tr = 70 kHz, P = 32 mW, w = 50 fim, 1/1.1,2 = 30 kHz 
(b) A = 790.06 ran, B = 4G, v rf = 3.295 MHz, ft r/ /27r = 270 kHz, P = 35 mW, w = 50 fun, 1/11,2 = 30 kHz 
(c) A = 800 nm, B = 4G, v r} = 2.78 MHz, fi r/ /27r = 25 kHz, P = 25 fiW, w = 1 fim, i/_li )2 = 30 kHz, Ax = 0.6w 



[14| | . The methods we describe may be used with or with- 
out additional tunnelling in the perpendicular directions 
(along spin independent potentials). Note that a low 
dimensional lattice can already have a larger U than a 
higher dimensional lattice due to tight confinement being 
possible along the perpendicular directions. 

In this work we confine our analysis to bosons, specifi- 
cally to 87 Rb. We choose to work in the F = 1 lower 
hyperfine state with three magnetic substates. With 
B r f and Bq perpendicular (Bq being the static magnetic 
field), there will be equal intensities of a + and o~ rf 
dressing fields. We assume a state independent (tt po- 
larised, or with large frequency detuning) optical lattice 
in the y and z directions, meaning that the difference in 
energy between magnetic substates is independent of y 
and z. We can then write, for these three substates, 

/y +1 (x) + <5 +1 n r/ /2 o \ 

H= n rf /2 V {x) Q rf /2 (2) 

\ o n P/ /2 v_i(x) + / 

with 5±i = AE±i_o/h±uj r f. The optical dipolc potentials 
V m p (r) can be shown to be [2(| 

7rc 2 r (l-Pgpm F 2 + Pg F m F \ 

VmAr) = { a Di + a D2 ) /(r) 

where Ap> 1 and A]j 2 are the frequency differences from 
the D lines, V = (f Dl + T D2 )/2, oj = (w Dl + 2w£> 2 )/3 
and P = ±1, for optical ct ± , tt polarisation. 

Using a Born-Oppenheimer-type approximation [2l| . 
the internal and external degrees of freedom of the atom 
can be decoupled as long as the kinetic energy of the 
atoms is much less than the energy spacing of the dressed 
levels. This is the case for typical fl r f, although Landau- 
Zener losses may occur as a result of this condition being 



weakly violated, as discussed below. Under this approxi- 
mation, the eigenvalues of the above matrix give the po- 
tential energy of the atoms in the corresponding dressed 
eigenstate. The rf coupling between the magnetic sub- 
states gives rise to avoided crossings, as shown in Fig. 

m 

We calculate the one dimensional Wannier functions 
for the lowest eigenstate of the Hamiltonian (Eqn. 
corresponding to the local ground state of the atoms. 
We diagonalise the Hamiltonian of a single potential pe- 
riod, and so find the eigenfunctions of the dressed poten- 
tial by using Bloch's theorem. We find the maximally 
localised Wannier functions at the lattice sites by recur- 
sively rephasing the eigenstates [22l j before summing over 
quasimomentum. From the lowest band Wannier func- 
tions we calculate the Bose-Hubbard parameters 

U = 2a s Ky/uJi_\u±_2 J \w {x)\ A dx (4) 

J=w^ E " eiqXr (5) 

q q 

with lo±i^2 the trapping frequencies in the perpendicu- 
lar directions, q the quasimomentum, E q the energy of 
the eigenstate with quasimomentum q, x r the separation 
between neighbouring lattice sites, and a s the scatter- 
ing length in the dressed state (for 87 Rb, approximately 
equal to the undressed scattering length as a s m a t [23j). 
The photon scattering rate is calculated from the Wan- 
nier functions and details of the dressed potential. 

We compare the on-site interaction energy between 
dressed and undressed potentials for a certain value of 
Hubbard J. We choose J/h = Jn m /h = 25 Hz so that 
the tunnelling time T tun = h/2zJ is Ttun — 10 ms (z 
being the coordination number). We choose this value 
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FIG. 2: (Colour online) Calculated properties of atoms in the ground state of the dressed potentials. The optical 
power P is the independent variable; dressing parameters v r $ and f2 r / are chosen subject to the conditions that the 
neighbouring wells are the same depth and that the tunnelling parameter J/h is 25 Hz. These conditions specify v r f 
and f2 r / for (b) and (c). For (a) the conditions only specify v r f 1 so Q r f can be chosen arbitrarily; two choices are 
shown. The dashed horizontal line gives the undressed U for comparison, a single value calculated for J/h = 25 Hz 
in the undressed er + or er~ potential. The scattering rates are from the dressed lattice beams only. Aside from P, 
i> r f and £l r f, the parameters of the potentials are as given in Fig. [1] 



as this is around the limit where slow interwell dynamics 
still occur on a typical experimental timescale; as J is de- 
creased further the phase transition will go from being an 
adiabatic to a non-adiabatic process. Interwell dynamics 
then cease, with the final purity of the Mott insulator 
depending on Uu m , the value of U when J = Ju m . 

The first potential we consider is given in Fig. Hal the 
'truncated sinusoidal' potential. This is an optical lat- 
tice at the 87 Rb, F = 1, mp = tune-out wavelength 
of 790.06nm [24J. The magnetic field magnitude and di- 
rection are chosen so that the mp = +1 state is detuned 
well above the other two states by the nonlinear Zecman 
effect. Experimentally plausible parameters are chosen. 
The lowest parts of the mp = — 1 potential intersect the 
flat mp — potential. Once dressed, the lowest adi- 
abatic potential differs greatly from a sinusoidal shape, 
with 'pockets' of strong confinement separated by an al- 
most flat potential. There is typically only a single bound 
state at each site. The enhancement of Uu m , for the cho- 
sen value of Jum, in this potential is given in Fig. [^a| 
for Q r f — 2tt x 70 kHz and 2ir x 10 kHz; the maximum 
increases we find are around 70% and 80% respectively. 
The enhancement factor is only weakly dependent on Q r f 
for a large range of rf power. 

The second potential we consider is given in Fig. Ilbl 
the 'interpolated sinusoidal' potential. Optically, it is 
very similar to the truncated sinusoid, but at a much 
lower magnetic field of a few Gauss, so the mp = ±1 
states are approximately symmetrical around the mp = 
state. In this case the rf dressing doubles the number 
of wells in the lattice. The enhancement of Uu m in this 
potential is given in Fig. l2b[ the typical increase is around 
150%. 

The last dressed potential (Fig. [Tc]) we consider is 
based on an array of independently addressable dipole 



traps[31|. We choose neighbouring dipole traps to have 
the same intensity but opposite a polarisation. Due to 
the rf dressing, the barrier between the potential min- 
ima can be considerable even when neighbouring Gaus- 
sian spots are optically unresolvable. The enhancement 
of Uum is shown in Fig. [2cl the maximum enhancement 
due to the dressing is around 80%. Note that, for poten- 
tials (b) and (c), the atomic spin adiabatically flips when 
a single atom tunnels from one site to its neighbour, en- 
abling sublattice addressability and readout. 

The results of these calculations show that the dressed 
Uii m is largely independent of the optical power, but 
is dependent on the dressing scheme. The results for 
the 'truncated' sinusoidal potential show that radically 
altering the shape of the periodic potential, to confine 
tighter while still allowing tunnelling, but without chang- 
ing the distance between wells, can only increase Uu m by 
a modest amount. However, Uu m becomes substantially 
larger if atoms are confined in the 'interpolated' sinu- 
soidal potential, with neighbouring wells separated by 
A/4. The interpolation method would therefore seem the 
most promising way to increase Uu m by rf dressing an 
optical lattice. 

The results show that the on-site interaction energy 
of atoms confined in a focussed Gaussian array can also 
be increased by the interpolation technique. The typi- 
cal values of Uu m in this case are not greatly less, for 
our chosen parameters, than for atoms in an undressed 
3D counterpropagating lattice. This raises the prospect 
of site-addressable highly number squeezed states in a 
hybrid focussed Gaussian / 2D optical lattice apparatus, 
with a comparable on-site interaction energy per atom to 
that observed in a Mott Insulator in a 3D optical lattice 
i. 

Even a modest increase in Uu m could be useful, 
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though, as the number of imperfections is likely to be 
a strong function of Uu m . Estimates of the thermal exci- 
tations would be nth ~ Ae~ (the J <^U,kT <^U limit 
of the Bose-Einstein distribution) and of non-adiabatic 

excitations n na w Be~ cu ; the latter from the form of 
the Landau-Zener avoided crossing transition probability 
in the limit J <C U » eoi, where eoi is the energy 
difference between the ground and first excited states of 
the many body system. 

We briefly consider how atoms may be loaded into 
the potentials (when J > U). Adiabatic loading of the 
'truncated' sinusoidal potential can be seen to be accom- 
plished by simply ramping the rf frequency, once atoms 
are trapped in the bare lattice. Adiabatic loading of 
the 'interpolated' sinusoidal potential and Gaussian ar- 
ray may be accomplished by ramping the laser intensity 
(with the rf already on), but the mechanism is less obvi- 
ous. In effect, during the final stages of the ramping of 
the lattice intensities, the atoms adiabatically delocalise 
between the two sets of potential minima as long as the 
rate of change of the offset between the two sublattices 
ej is much less than the other Hamiltonian parameters. 

It is important to note that spin-changing collisions 
and other 2-body collisional loss rates are greatly s up- 
pressed for atoms trapped in the lowest eigenstate [ll|, 
[2(1 \2l\ . This is because a collision which changes the 
internal states of these atoms needs an input of energy 
which is significantly greater than the energy available 
from kinetic or potential energy (of order J or U) . The 
other loss mechanism is Landau-Zener losses arising from 
atomic motion at the avoided crossings, which have dom- 
inated the only rf dressed optical potential experiment 
carried out to date [id ]. However, Landau-Zener losses 
become negligible when the atoms are trapped in the low- 
est eigenstate, as the transition from the lower to the up- 
per dressed level needs energy which is required to come 
from the kinetic energy of the atoms; this kinetic energy 
~ J is typically very much less than f2 r /, so Landau-Zener 
losses can be expected to be negligible for atoms in the 
lowest trapped state [281 ]. 

A remaining loss mechanism, specific to the truncated 
lattice, is the coupling of weakly bound states to the con- 
tinuum. This can only occur when the highest energy 
state in the lowest band becomes greater than the bind- 
ing energy of the site. When mU and J are significantly 
less than the binding energy of the site, as is the case 
in the examples given above (assuming around one atom 
per site), this loss mechanism is, to a large extent, sup- 
pressed; however with higher numbers of atoms per site 



there is a chance that some atoms will be ejected from 
the lattice. 

The calculations shown do not specify whether 
tunnelling is proceeding in the perpendicular (spin- 
independent) directions; J along the dressed lattice di- 
rection is independent of J in the other directions. The 
lattice parameters in the other directions only influence 
U through a multiplicative factor in Eqn. 2] Therefore, 
when there is tunnelling in all directions in such a lattice, 
the enhancement factor for U arising from dressing along 
a single direction will be the same whether or not there 
is also tunnelling in the perpendicular directions. Thus 
this technique may be used in a 3D lattice (composed of 
ID spin-dependent and 2D spin-independent lattices), to 
enhance Uu m by the same factor as in purely ID lattices. 

If spin-dependent potentials are present along more 
than one direction, the 2D or 3D rf dressed potential 
cannot be simply expressed as a sum of ID potentials; 
one consequence is that the potential minima may occur 
along lines or surfaces rather than at points. A short 
discussion may be found in Ref. [14| ; we do not consider 
these cases further in this paper. 

A practical complication with this scheme, in common 
with other rf dressed optical potential schemes, lies in 
the use of magnetic field sensitive transitions. The typ- 
ical magnetic field drift in a laboratory environment is 
of order ^lmG [29(; this would manifest as a sublattice- 
dependent energy offset of order ^lkHz. Clearly we need 
this site-dependent offset to be less than U; practically 
this would mean using a magnetic shielding technique 
such as mu-metal cladding, which can decrease these am- 
bient fields by a factor of around 100 [301 ] . 

In conclusion, we show that substantial enhancement 
of on-site interactions may be achieved by rf dressing op- 
tical potentials. We study three cases of interest, and 
show that it is possible to enhance nonlinear parameters 
by more than a factor of 2, with negligible Landau-Zener 
losses, which should greatly improve the purity of the 
resulting Mott insulator state. We find that decreasing 
the distance between neighbouring wells has a greater ef- 
fect on the limiting on-site interaction than modifying the 
form of the potential between existing wells. In summary, 
the techniques described have the potential to make com- 
plex quantum states with neutral atoms purer and faster. 
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